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Abstract: 

We reconsider the prediction of the semileptonic branching ratio of B mesons, using the 
recent calculation [|l| of the radiative corrections with account for finite quark masses 
in nonleptonic decays and taking into account I /ml corrections. For the semileptonic 
branching ratio we obtain Bsl = (11.8±1.6)% using pole quark masses and Bsl = (11.0 ± 
1.9)% using running MS quark masses. The uncertainty is dominated by unknown higher 
order perturbative corrections. We conclude that the present accuracy of the theoretical 
analysis does not allow to state a significant disagreement with the experimental results. 
However, our re-analysis of the decay b —>■ cos yields an increase of (35 ±11)% due to next- 
to-leading order corrections including mass dependent terms, which further emphasizes the 
problem of the average charm quark content of the final states in B decays. 
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1. The theoretical description of inclusive weak decays of heavy hadrons has made 
considerable progress over the recent years, see for a review. It could be shown that in 
the limit of infinite heavy quarks the decay rate coincides with that of the corresponding 
free quark decay; corrections to that limit are of nonperturbative origin and suppressed by 
at least two powers in the heavy quark mass 0. Today there is increasing confidence that 
QCD predictions of heavy particle decays rest on a firm theoretical foundation. In view 
of these apparent advances and with the availability of new and more precise data on the 
semileptonic branching ratio Bsl of the B meson [Q, the long felt discrepancy between its 
measured and its predicted value becomes more and more baffling P] . Over the last years, 
the main efforts were concentrated on the determination of the size of nonperturbative 
power-suppressed corrections to the free quark decay, which, however, turned out to be 
small, of natural size ~ IGeV^/m^ ~ 5%, and cannot explain the experimental value of 
Bsl- Thus, it seems timely to place more emphasis on perturbative radiative corrections 
to the free quark decay, which since the well-known analysis by Altarelli and Petrarca 
have not been receiving proper attention in the literature. In this letter we update 
the theoretical prediction of the semileptonic branching ratio of B mesons using a recent 
calculation of the charm quark mass dependence of radiative corrections to nonleptonic 
decays. In addition we re-analyze the decay rate r{b —>■ ccs), taking into account the 
quark mass dependence of radiative corrections and the contributions of penguins. The 
semileptonic branching ratio is then evaluated using both pole masses and running quark 
masses. The latter procedure was advocated in 0] on the evidence of the cancellation of 
renormalon singularities 0. Finally, we discuss shortly the problem of fixing a proper 
renormalization scale in heavy quark decays. 

2. The semileptonic branching ratio of B mesons is defined by 

o ^ nS - Xev) 

T.^=e,^,,r^{B Xiug) + r{B X^) + r{B + r(rare decays) ' ^ ' 

The heavy quark expansion (HQE) allows to relate the inclusive decay rate of a B meson 
to that of the underlying b quark decay process, apart from 1 / corrections: 

T{B ^X) = T{b ^x) + Oil/ml). (2) 

The power-suppressed correction terms to the total inclusive widths of both semi- and 
nonleptonic decays were calculated in 0, |§] . 

For the free quark decay rates we introduce the following notations: 

T{b ^ ciu) = To PH(xe, xe, 0) /(x^, x,, 0), (3) 
T(b cud + cus) = 3Tq PII(xc, 0, 0) //(/i) J(xc, /i), (4) 

r(6 — i> CCS + ccd) = 3Tq PH(xc, Xc, Xg) i^{xc, Xs, /i) K{xc, Xs, /i). (5) 

Here Fq is defined by Fq = G'jp\Vcb\'^ ml / {192n^) . PH(xi, 0:2, X3) is the tree level phase space 
factor of the decay b ^ qi + W qi + q2 + Qs] for arbitrary masses Xi = nii/mb it is given 
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by: 



/ds 
— {s-xl- xl) {l + xl~ s) w{s, xl, xl) w{s, xj, 1) (6) 



ds 

{X2+X3) 

with 



w{a, b, c) = (a^ + + - 2ab - 2ac - 2bcy^\ (7) 

The functions rj and k, contain the leading-order QCD corrections to the nonleptonic rates 
b —>■ cuq and b — >• ccq, respectively. In particular, r] is given by |T0| 



with /5o = 11 — 2nj/3 for Uf running flavours, = 5 in our case. The expression for K{fi) 
is given below. Finally, /, J and K contain the next-to-leading QCD corrections to the 
decay rates. The function I can be written as 

2 as 

I{xi,X2,X3) = 1 + - —g{Xi,X2,X3), (9) 
6 TT 



where g has been calculated in Ref. W^ for arbitrary arguments in terms of an one- 
dimensional integral. Analytic expressions are available for the special cases g{xi, 0, 0) |[T2| 
and g{0,X2,0) [Q]. The complete calculation of J{xc,fi) was flrst done in [T|, while J(0,/i) 



is also available from [|T^; some of the terms for arbitrary Xc have been also calculated in 



| rT| . The function K is not yet known completely; we will discuss it below. 

Summarizing existing calculations of the radiative corrections, we give the numerical 
values of g{xc,0,0), g{xc,Xr,0) and J{xc,mh) in Table 1. The numbers are evaluated for 
«s(^z) = 0.117, i.e. = 312 MeV, and at the renormalization scale n = rrib = 4.8 GeV. 
With these parameters, we find i]{mi,) = 1.10. 

3. As explained above, all the decay rates entering the semileptonic branching ratio (0) 
are known to next-to-leading order in the strong interaction including final state particle 
mass effects, except for T{b ^ ccs + ccd) and the rare decays. Whereas the latter can 
safely be neglected, the channel b —>■ ccs deserves a closer consideration. In addition to the 
contributions studied in Ref. [|l|, where a 30% increase of the decay rate b —>■ ccs by radiative 
corrections was obtained^, we take into account the dependence of these corrections on the 
s quark mass and discuss the contributions of penguins. 

The leading order decay rate can be written as 



r{b^ccs)\^^ = 3ro|K«|'PH(xe,x„a;,)<j5]c2(/i) + 2 



^Ci{fl)c2{fl) + Ci{fl)c3{fl) 



^See also 1 14 
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+ ^ Clifi)c4{fl) + ^ C2(Ai)c3(/i) + C2(Ai)c4(/i) + ^ C3(/i)c4(/i) + ^ C5(/i)c6(/i) 



Ci(^)c5(^) + ^ Ci(/i)c6(/i) + ^ C2(/i)c5(/i) + C2(/i)c6(/i) 



+C3(/i)c5(/i) + ^ C3(/i)c6(/i) + ^ C4(Ai)c5(/i) + C4(Ai)c6(Ai) 



(10) 



= 3To\Vcs\'^PB.{Xc,Xc,Xs) K{Xc,Xs,fi). (11) 

The coefficients Ci{fi), 1 < ^ < 6, are the leading order Wilson-coefficients multiplying 
the operators Qi in the effective Lagrangian and can be found in tabulated form in [p!5| . 
Unlike the expression in Ref. 0], our Eq. (|TUp also takes into account the interference of 
four-quark operators having the usual (V — A)(S)iV — A) structure with penguin operators 
of the structure {V — A) ® (V + A). These interference terms are explicitly of order x^ and 
enter the decay rate with a weight-function /, given by 

1 

/(xc,Xc,x,) = — — / ds^w{s,xl,x^^)w{l,s,xl){s + x^,-xl){l + s-xl). 

(12) 

For reasonable quark masses Xc = 0.3 and Xg = 0.04 we find / = 0.24 and K,{fi = rrih = 
4.8 GeV) = 1.07. Neglecting the penguin-contributions, i.e. for Cj(yu) = for i > 3, ^{mi,) 
coincides with ri{m}y) = 1.10, so that the penguins interfere destructively and reduce the 
decay rate by ~ 3% similarly to what was observed in 0. 

In next-to-leading order, the decay rate can be written as in Eq. (|^), where K is defined 

by 

6 

k[Xc, Xs, fJ')K{Xc, Xg, fi) =: fiji^Xc, Xs)Ci{/J,)Cj{fi)dij(^Xc, Xg, fi) , (13) 

the weight-factors fij being given in (0), whereas the dij have the structure 
d., = l + h,^ + r., ^s{m^) - ^sif^) ^ , 

71 TT 

The terms rjj contain matching-coefficients and two-loop anomalous dimensions of the 
operators Qi and can be obtained from [15|. The terms ku and k22 can be obtained from 
in, |lT], including all dependence on Xc and Xg, likewise for Xc = Xg = and, partly, also 
in dependence on Xc and Xg. The other terms are not known. Nevertheless, the knowledge 
of these three coefficients allows a rather accurate determination of the decay rate to next- 
to- leading accuracy: for the unknown d^, we most conservatively assume < d^ < 2, 
which corresponds to \kij\ < 15 for Og = 0.2. For di2, we replace the uncalculated term, 
He{Xc,Xg) in the notation of [|I|, by its corresponding value for only one massive c quark, 
Ge{xc). We estimate the error introduced by this procedure by AH^, ^ 2\Ge{xc) — ^6(0)1. 
The values of the relevant kij are given in Table 2, together with the functions n and 
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the latter one yielding the increase of r(6 ccs) due to next-to- leading order QCD 
corrections. 

At this point, it is worthwhile to emphasize that the dependence of the decay rate on 
the s quark mass is rather weak. While the phase space factor is considerably reduced by 
including the strange mass 0], this effect turns out to be to a large extent compensated 
by the increase of radiative corrections. For Xc = 0.3 we find that a strange quark mass 
rris = 200 MeV, Xg = 0.04, reduces the decay rate by 1.5% only, which is smaller than the 
effect of the penguins. 

Taking everything together, we find that for Xc = 0.3 and Xg = 0.04, next-to-leading 
order radiative corrections increase r(6 — ccs) by (35 ± 7^j)%, where the first error is a 
very conservative estimate of the unknown parts of the next-to-leading corrections and the 
second error comes from a variation of the renormalization scale /z within mh/2 < /i < 2mfe. 

4. We have now all ingredients at hand to evaluate Bsl- Before doing so, however, 
let us make some remarks about the nonperturbative corrections entering the decay rates. 
They can be expressed in terms of two hadronic matrix elements, Ai and A2. Whereas A2 
is directly related to the observable spectrum of beautiful mesons, 

A2^^(m|. -m|) = 0.12GeV2, (15) 

the quantity — Ai/(2mfc), which can be interpreted as the average kinetic energy of the b 
quark inside the B meson, is only difficult to measure. At present, only a QCD sum rule 
calculation is available, according to which Ai = —(0.5 ± 0.1) GeV^ [l^. For a summary 



of the discussion about Ai we refer to JT^- The formulas for the decay rates including 
power-suppressed corrections are given in 

We next have to fix the quark masses that enter the decay rates. For the strange quark 
mass we use = 0.2 GeV; as emphasized in the last section, the decay rates are not very 
sensitive to this parameter. As for the heavy quark masses, we make use of the fact that 
in the framework of HQE the difference between and is fixed: 

Ai + 3A2 / 1 1 \ I \ 

mb-mc = mB-mD^ + C ^ . (16) 

2 \mh rricJ \'^qJ 

For TTifo we use = (4.8 ± 0.2) GeV. Varying the renormalization scale within the range 
mb/2 < /i < 2mfc, we find 

Bsl = (11.8 ± 0.8 ± 0.5 ± 0.2 ± 0.21?;^)% , (17) 

which is our main result. Here the first error comes from the uncertainty in rrib, the second 
one from the one in as{mz) = 0.117 ± 0.007, the third one from the uncertainty in Ai and 
the fourth from the uncertainty in r(6 ^ ccs). The last error comes from the variation of 
the renormalization scaleFl. 



^ Larger values of mh and (or) of the normalization scale generally yield a larger Bsl , while the increase 
of as tends to lower the branching ratio. 
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The error stemming from the uncertainty in n is rather big and shows that higher order 
perturbative corrections are important. We thus feel motivated to evaluate Bsl also in 
a different scheme, using running short-distance masses, e.g. MS masses. This procedure 
has also been advocated in connection with the cancellation of renormalon contributions 
0]. In order to keep the formulas scheme-independent at 0{as), the phase-space has to be 
modified according to 

ml PH(a;„ 0, 0) — ml PH(a:„ 0, 0) |l + ^ f f - 5 In ^ - In d\nFRjx^,0,Oy 

where Xi denotes running quantities evaluated at the scale fi. With this substitution, we 
obtain 

Bsl = (11-0 ± 0.6 ± 0.8 ± 0.2 ± 0.2+^;°)% (19) 

In Table 3 we give a comparison of theoretical predictions for Bsl using different 
approximations. The main result of our analysis is that the prediction of Altarelli and 
Petrarca is lowered by more than 1%. It is clearly visible that the main effect comes from 
taking into account the quark mass dependence of radiative corrections calculated in [Q, 
while the nonperturbative l/mf corrections result in a ~ 0.2% decrease, in agreement with 

§■ 

The results shown above clearly demonstrate that the main theoretical uncertainty in 
the predictions for B decays comes from the scale and scheme dependence of the results. 
As always, the only consistent way to reduce both is to make a next-to-next-to-leading 
(NNLO) calculation of the decay rate including Ola"^) corrections, which is a formidable 
enterprise. Lacking this calculation, one is bound to make rather crude estimates of the 
possible higher-order corrections using some particular prescription to fix the scale. Among 
these, the BLM prescription seems to us to be the only one that is physically motivated. 



We now discuss in short its possible outcome on the B decay widths. 

The idea underlying the BLM approach is that the major part of higher-order radiative 
corrections originates from the necessity to evaluate Feynman diagrams with the running 
coupling at the scale of the gluon virtuality and can be traced by a relatively simple calcu- 
lation of the diagrams with an extra fermion bubble in the gluon line. The corresponding 
calculations have been done recently [|l^, ^ and indicate that the natural scale in the 
radiative corrections in B decays is significantly smaller than m?,. In particular, neglect- 
ing the c quark mass, Ref. ||20| finds hblm = O.OTm^ for the radiative corrections to the 



semileptonic decayj^ while for the final state interaction of quarks in the nonleptonic de- 
cays fJ,BLM ~ 0.32mb is found, as indicated by the studies of the r lepton hadronic decay 
width. In fact, the particular scale entering the radiative corrections to the semileptonic 
width turns out to be not very important for the problem of the semileptonic branching 
ratio, since these corrections cancel to a large extent between numerator and denominator 

Note that this result is contrary to expectations in [|j , where the choice of a low scale in the radiative 
corrections was criticized. 
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of Eq. (|1|).Q The main problem in applying the BLM scheme to nonleptonic B decays is 
that the radiative corrections in NNLO do not factorize into "semileptonic" and "final- 
state-interaction" parts. On may thus suspect that a large part of the radiative corrections 
comes from other types of diagrams than those considered in the BLM method (see the 



discussion in |T^). In addition, it is not clear how to apply this approach consistently to 
processes where the relevant operators possess a nontrivial anomalous dimension. Still, 
we believe that the low scales indicated by the BLM prescription are more natural in B 
decays. Thus the choice /i = 2mb adopted above as one extreme case is in fact rather 
unlikely, while /i = mb/2 is presumably more relevant. Adopting this scale, our result for 
the semileptonic branching ratio in (0) becomes Bsl = (10.5 ± 1.4)%, in perfect agree- 
ment with the experimental value = (10.4 ± 0.4)% [^]. Summarising, we conclude 
that there is no evidence for any disagreement between the experimental data and the 
theoretical prediction for the semileptonic branching ratio of B mesons. 

The situation is not so clear, however, with the charm content in the final states. With 
the 35% increase of the b ccs rate induced by taking into account the c quark mass in 
the radiative corrections, this problem is strengthened. From our analysis we get 

(ric) = 1.28 ±0.08, (20) 



which is to be confronted with the experimental result (ric)'^^^ = 1.04 ± 0.07 [^. We are 
not aware of any natural theoretical possibility to lower the value given in (|20|), unless the 
c quark mass is much larger than expected, which would confiict, however, with the heavy 
quark expansion of the meson masses, Eq. (0). 



^In the case at hand, where the BLM prescription indicates a very low scale, we find it more appropriate 
not to change the scale fi = mf, to fi ~ O.OTm;,, but rather to use the calculation in [ pO| as an explicit 
estimate of the a^lmi,) correction. 
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Tables 





n(r 01 




1 i T TYl L 1 
V C' '"'Oj 


n 

yj 




— 37 


1 OOQ 


n 1 




—9 8Q 


1 096 


0.2 


-2.84 


-2.42 


1.046 


0.3 


-2.51 


-2.08 


1.063 


0.4 


-2.23 


-1.81 


1.077 


0.5 


-2.01 


-1.61 


1.088 


0.6 


-1.83 


-1.45 


1.097 


0.7 


-1.70 




1.105 


0.8 


-1.59 




1.113 


0.9 


-1.53 




1.123 


1 


-1.50 







Table 1: Next-to-leading order corrections to the semileptonic b quark decay rates and 
the decay h — > cud as functions of = mc/mi,. Parameters: iJ, = mi, = 4.8 GeV, = 
312 MeV, corresponding to as{mz) — 0.117; x^- — mr/mij. Note that J(l,/^) diverges like 
~ ln(l - Xc). 





K,{X(., Xsi 


kii 


^'12 (/U = rrih) 


^^22 


K{xc,Xs,mh) 





1.054 


-1.33 


-7.59 ±0.01 


-1.26 


1.02 ±0.05 


0.1 


1.056 


-0.05 


-6.65 ±0.07 


-0.35 


1.09 ±0.06 


0.2 


1.062 


2.53 


-4.97 ±0.20 


1.23 


1.20 ±0.06 


0.3 


1.069 


6.69 


-2.64 ±0.57 


3.41 


1.35 ±0.07 


0.4 


1.077 


15.68 


1.24 ±0.96 


7.09 


1.62 ±0.09 



Table 2: The leading and next-to-leading order corrections to the nonleptonic decay b — > 
CCS. The errors rely on a conservative estimate of the unknown parts of the next-to-leading 
order terms, mostly due to penguin contributions. The last column gives the increase of 
the decay rate r(6 — > ccs) in next-to- leading order including finite c and s quark effects in 
the radiative corrections. The input parameters are the same as in Table 1; Xg — 0.04. 
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Parton Model 


HQE § 


HQE [this work] 




pole masses 


pole masses 


pole masses MS masses 


0.110 


0.133 


0.132 


0.123 0.117 


0.117 


0.130 


0.128 


0.118 0.110 


0.124 


0.125 


0.123 


0.113 0.102 



Table 3: Theoretical predictions for the semileptonic branching ratio Bsl depending on 
as(^z)- Input parameters: rrib = 4.8 GeV, rric = 1.33 GeV (pole masses), corresponding 
to Ai = —0.5 GeV^, = 0.2 GeV. Renormalization scale: /x = rrib. 
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Erratum: 

Theoretical Update of the Semileptonic 

Branching Ratio of B Mesons 
[Phys. Lett. B 342 (1995) 362] 

E. Bagan, Patricia Ball, V.M. Braun and P. Gosdzinsky 

In the third line in Eq. (10) on page 364 there is a sign error: —2f(xc,Xc,Xs) should 
read +2f{xc,Xc,Xs)- In addition, we have found an error in the computer program, which 
affected the average charm content Uc and the scale dependence of the results in the MS 
scheme. We take this opportunity to incorporate the complete results for the quark mass 
dependence of the radiative corrections to 6 ^ ccs calculated in [22]. The corresponding 
update of our Table 2 on page 365 is given in Table 2 in [22] . 

The numerical impact of these corrections is marginal: Eqs. (17) and (19) on page 366 
should read: 

BsL = (12.0 ± 0.7 ± 0.5 ± 0.2l°|)% , (17) 

BsL = (11.3 ± 0.6 ± 0.7 ± 0.21°:?)% . (19) 

Table 3 on page 367 has to be replaced by the Table given below. 

Since the problem of the average charm content is receiving increasing attention (see, 
e.g. [23]), we give the corrected result for ric in a somewhat expanded form. Eq. (20) on 
page 367 is to be substituted by 

nc = 1.24 ±0.05 ±0.01, (20) 

which shows the result in the OS scheme. The first error comes from the uncertainty in 
rrib — (4.8 ± 0.2) GeV, the second one from the uncertainties in the remaining parameters. 
In the MS scheme we get 

ric = 1.30 ±0.03 ±0.03 ±0.01, (20') 

where again the first error comes from the uncertainty in the quark masses, the second one 
is due to the variation of as, and the third one comprises the remaining uncertainties. 

We have added a figure showing the charm content versus the semileptonic branching 
ratio, cf . [23] , obtained by relaxing the constraint on the quark masses following from the 
heavy quark expansion in Eq. (16) on page 366 and allowing for a larger range of the ratio 
rric/mi,. Note that rric/mb is scale-independent; both ric and Bsl are functions of rric/mb, 
/J, and Q!s(/x). 

Acknowledgement: We thank G.Buchalla and M. Neubert for pointing out the errors. 
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Parton Model [7] 


HQE [5] 


HQE [this work] 


asimz) 
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pole masses 


pole masses MS masses 


0.110 
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0.132 


0.124 0.120 


0.117 


0.130 


0.128 


0.120 0.113 


0.124 


0.125 


0.123 


0.114 0.105 



Table 3: Theoretical predictions for the semileptonic branching ratio Bsl as a function 
of as{mz). Input parameters: rUb = 4.8 GeV, rric = 1.33 GeV (pole masses), corresponding 
to Al = —0.5 GeV^, nig = 0.2 GeV. Renormalization scale: fj, = mi,. 
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0.1 0.11 0.12 0.13 0.14 

BsL 

Figure 1: The charm content ric vs. Bsl- Sohd hnes: theoretical predictions in the OS 
scheme for 0.23 < rric/mi) < 0.33, dashed hnes: the same in the MS scheme for 0.18 < 
'nT-ci'i^c) /i^b{i^b) < 0.28. Shaded area: theoretical predictions in the OS scheme with rric 
obtained from Eq. (16) and varying TJib, Ai and fi within the range of values given in the 
text. The experimental data point is taken from [24]. 
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